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Discretized landscapes can be mapped onto ranked surfaces, where every element (site or bond) 
has a unique rank associated with its corresponding relative height. By sequentially allocating these 
elements according to their ranks and systematically preventing the occupation of bridges, namely 
elements that, if occupied, would provide global connectivity, we disclose that bridges hide a new 
tricritical point at an occupation fraction p = pc, where pc is the percolation threshold of random 
percolation. For any value of p in the interval pc < P < 1, our results show that the set of bridges 
has a fractal dimension dsB ~ 1.22 in two dimensions. In the limit p ^ 1, a self-similar fracture 
is revealed as a singly connected line that divides the system in two domains. We then unveil how 
several seemingly unrelated physical models tumble into the same universality class and also present 
results for higher dimensions. 



Any real landscape can be duly coarse-grained and rep- 
resented as a two-dimensional discretized map of regular 
cells (e.g., a square lattice of sites or bonds) to which 
average heights can be associated. This process is exem- 
plarily shown in Figs. [TJa)-(c). As such, the concept of 
discretized maps has been considered as a way to delimit 
spatial boundaries in a wide range of seemingly unre- 
lated problems, ranging from tracing water basins and 
river networks in landscapes [T]-[5] to the identification 
of cancerous cells in human tissues [G', T , and the study 
of spatial competition in multispecies ecosystems [8[ |9] . 
Moreover, previous studies have shown that cracks or 
surviving paths through discretized maps possess a uni- 
versal fractal dimension which can be physically realized 
in terms of optimal paths under strong disorder [10 13 , 
optimal path cracks JHITl], loopless percolation [16 , 17 , 
or minimum spanning trees p!8H23] . Here we show that 
all these problems can be understood in terms of the same 
universal concept of fracturing a ranked surface. 

We start by defining a ranked surface. Given a two- 
dimensional discretized map of size L x we generate a 
list containing the heights of its elements (sites or bonds) 
in crescent order, and then replace the numerical values 
in the original map by their corresponding ranks. As de- 
picted in Fig. [ijd), the result is a ranked surface. The 
process of fracture generation is rather simple. Once the 
ranked surface is obtained, we sequentially occupy the 
elements of an empty lattice with the same size follow- 
ing the crescent rank order of the corresponding elements 
(i.e., in the same position) on the ranked surface. During 
each step of the allocation process, only bridges, identi- 
fied as those lattice elements which, once occupied, would 
create a spanning cluster (i.e., a globally connecting clus- 
ter) [24], are never occupied. These elements will even- 
tually form a macroscopic fracture. 

In Fig. |2] we show the evolution of the fracture line 
on a large ranked surface with the fraction of occupied 
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bonds p. As displayed, the lattice is initially seeded by 
a set of disconnected bridge elements at low values of p, 
while for p ^ 1 the fracture finally emerges towards a 
singly connected line that divides the system in two. As 
we show later in this article, our results reveal that this 
line is fractal with dimension (Ibb ^ 1.22. Interestingly, 
this value is statistically identical to the dimensions of 
fractures generated from different models previously in- 
vestigated p!QHT3] . However, at the percolation threshold 
value of the classical random percolation model [24l [25] , 
p = Pci the set of bridge bonds appearing in any con- 
figuration of our model should be identical to the set of 
the so-called anti-red bonds in random percolation |26] . 
As first proposed by Coniglio [26 and numerically veri- 
fied by Scholder [27], at p = this set is also fractal, 
but with dimension 1/v (= 3/4 in 2D), where v is the 
correlation length exponent. 

Motivated by this substantial change in the fractal be- 
havior at distinct stages of our fracturing process, in what 
follows we address how the set of bridges scales with the 
fraction of occupied bonds and the system size at p = 
analogously to a theta point [28^,30 , while for all values 
of p above it has a fractal dimension dBB- More- 
over, we introduce a new tricritical crossover exponent, 
which we study up to dimension six, the upper-critical 
dimension of percolation. 

I. RESULTS 

We performed simulations of fracturing ranked surfaces 
on square lattices. It is worth noting that, despite the 
similarities with random percolation, the suppressing of 
connectivity poses a statistically different problem. For 
example, while for p ^ 1 there is only a single config- 
uration in random percolation (all bonds occupied), in 
ranked percolation there are A/"!, evenly weighted, possi- 
ble configurations, where N is the total number of bonds. 
In classical percolation the total number of configurations 
is 2^. Figure [3] shows the dependence of the number of 
bridge bonds Nbb on system size, for different fractions 
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FIG. 1. The generation process of a ranked surface. The landscape in (a) is coarse-grained to the low-resolution system of 
8x8 shown in (b), and then represented as a discretized map of local heights, as depicted in (c). By ranking these heights in 
crescent order, one obtains the ranked surface in (d). In fact, the landscape shown in (a) is a high resolution synthetic map 
obtained from a fractional Brownian motion simulation based on the Fourier filtering method [151 1411446] . 



of occupied bonds, namely, p = Pc = 0.5, p = 0.51, and 
p = 0.8. As expected, at the percolation threshold of 
classical percolation {p = p^)^ the number of bridge bonds 
diverges with system size as Nbb ^ L^l^ ^ where v is the 
correlation length exponent, with v = 4/3 in 2D; while 
for p = 0.8, Nbb ^ L^"""" ^ with cIbb = 1.215 ± 0.003. 
The latter is in fact observed at any p > pc. We found 
the same result (see Appendix) for site percolation and on 
other lattices (star, triangular, and honeycomb), which 
provides strong evidence for the universality of this ex- 
ponent. For p ^ Pc, like p = 0.51 we can observe, as 
depicted in Fig. |3j a crossover between the two differ- 
ent regimes. The inset of Fig. |4] shows Nbb, rescaled 
by L^^^, as a function of p, for different system sizes. 
The number of bridge bonds grows with p, such that, 
Nbb ^ {p- Pc)^, where C = 0.50 ± 0.03 is a novel expo- 



nent, which we call bridge-growth exponent. The overlap 
of the different curves confirms that the fractal dimension 
of the bridge bonds above Pc is (Ibb, for all p > Pc. This 
result differs from classical percolation where fractality 
is solely observed at criticality j24l [25] while, above Pc^ 
bridge bonds are only observed for finite systems [26^ . 



For polymer chains, at high temperatures, the ex- 
cluded volume prevails over attractive forces and the 
chain can be described as a self-avoiding walk. When 
the temperature is reduced, the attractive forces become 
relevant leading, at a theta-temperature, to a new ex- 
ponent at the crossover between two dimensions j28H3Q] . 
Analogously, in ranked percolation the fractal dimension 
of the bridge bonds is ai p = Pc, between (Ibb above 
Pc and zero below Pc. For the tricritical scaling we verify 



FIG. 2. Snapshots of the fractal set of bridges in two (hne) and three (surface) dimensions. For 2D four stages are seen (from 
left to right): p = Pc (black), p = 1.0 Ipc (blue), p = 1.05pc (green), and p = 1 (red), while for 3D only the final set is shown. 
We considered in 2D a square lattice with 1024^ sites and, in 3D, a simple-cubic lattice with 512^ sites. The fractal dimension 
is dBB = 1.215 =b 0.003, in 2D, and dsB = 2.50 ± 0.02, in 3D. 
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FIG. 3. Crossover in the size dependence. Number of bridge 
bonds, Nbb, for different p, namely, p — Pc — 0.5 (circles), 
p — 0.51 (stars), and p = 0.8 (triangles). The solid lines stand 
for the best fit. At p — pc, as conjectured by Coniglio [26] , 
Nbb ^ L^/^, where u = 4/3 is the critical exponent of the 
correlation length in 2D. For p > Pc, the number of bridge 
bonds scales with L^^^ . A crossover between the two regimes 
in system size is observed (stars) for p in the neighborhood of 
Pc. Systems of size have been considered, with L ranging 
from 32 to 4096. All results have been averaged over 10^ 
samples. Error bars are smaller than the symbol size. 



the following ansatz, 

Nbb = l'^'':f[{p-Pc)l'], (1) 

where T[x] ^ for x 7^ 0, and is nonzero at x = 0; and 



FIG. 4. Tricritical scaling, crossover, and data collapse. 
Number of bridge bonds, Nbb, as a function of the frac- 
tion of occupied bonds, p, for 2D with different system 
sizes L = {256,512,1024,2048,4096}. The scaling function 
given by equation ([T]) is applied, with = 0.93, obtaining 
C = 0.50 ± 0.03. In the inset, Nbb has been rescaled by 
L^BB ^ where (Ibb — 1.215. All results have been averaged 
over 10^ samples. 

the exponent is the crossover exponent. Therefore, the 
following relation is obtained, 

e = c' (d-^-^y (2) 

where (p = d — dsB- In the main plot of Fig. |4] we see, 
for 2D, the scaling given by equation ([T]), with = 0.93. 
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TABLE I. Table of exponents. Values of the exponents C 
and Lp up to dimension six. With increasing dimension, the 
C exponent converges to ^ — 1.5±0.7 and Lp goes to zero, 
revealing that the set of bridge bonds is dense. 



d 


c 




2 0.50 ± 0.03 0.785 ± 0.003 


3 


1.0 it 0.1 


0.50 ±0.02 


4 


1.3 ±0.5 


0.26 ±0.08 


5 


1.4 ±0.6 


0.1 ±0.2 


6 


1.5 ±0.7 


0.0 ±0.1 



The results above disclose a tricritical Pc below which 
the fraction of bridges in the bridge line vanishes in the 
thermodynamic limit. This is identical to a minimum 
height in the bridge line, H = /imin, in the context of 
landscapes (Figs. [TJa)-(b)). For a cumulative distribu- 
tion function of heights P{H < h)^ the minimum is 
given by P{H < /imin) = Pc Note that, for uncorre- 
lated landscapes, regardless the distribution of heights, 
the set of bridges only depends on their position in the 
rank. To observe the new set of exponents on discretized 
landscapes (Fig. [T](c)), Nbb is the number of sites in the 
bridge line with both neighbors (one at each side) having 
height lower than /i, where P{H < h) = p. 

To study the dependence of exponents ( and (p on the 
spatial dimension, we analyze the same problem up to 
dimension six. On a simple-cubic lattice {3D)^ above 
Pc^ the set of all bridge bonds has a fractal dimension 
dsB = 2.50 ± 0.05 and grows with C = 1-0 ± 0.1 (see Ap- 
pendix). Figure [s] shows the size dependence of Nbb^ in 
the limit p = 1, for lattices with size L^, where 2 < d < 6 
is the spatial dimension. In the inset, we plot the expo- 
nent as a function of d. Since the set of bridge bonds 
blocks connectivity from one side to the other, its fractal 
dimension must follow d — 1 < dsB ^ i.e., < (f < 1. 
With increasing dimension cp decreases. At the upper- 
critical dimension of percolation, dc = 6, (f = 0.0 ±0.1 
and the set of bridge bonds becomes dense having the 
spatial dimension d. Table |l] summarizes the exponents 
for dimensions 2 to 6. The bridge-growth exponent grows 
with d and converges toC = l-5±0.7at the upper-critical 
dimension. For (i > 6, above the critical dimension of per- 
colation, the exponent cp remains zero and the dimension 
of the set of bridges is equal to d. This can be understood 
from the fact that above dc one has an infinity of span- 
ning clusters [31 and thus many more possible bridges. 
Since the dimension of the set of bridges ai p = Pc is l/u 
and above Pc is the crossover exponent increases with 
d^ with the relation given by equation ([2|, where cp = 0. 

Since one can interchange occupied and empty bonds, 
there exists a symmetry between bridges and cutting 
bonds (red bonds), so one can raise the question of what 
happens when bonds are removed from a percolating sys- 
tem with the constraint that connectivity cannot be bro- 
ken. Initially all bonds are occupied and at the end, 
since cutting bonds are never removed, a line of cutting 
bonds is obtained, which we denote here as the cutting 




FIG. 5. Dependence on the spatial dimension and approach 
to the mean- field limit. Size dependence of the number of 
bridge bonds, Nbb, in the limit p = 1, up to dimension six. 
The solid lines stand for the best fit. The mass of the shortest 
path in loopless percolation (SP) in 2D is also included for 
comparison. Results have been averaged over lO'^ samples for 
2D and 3D, and 10^ samples for higher dimensions. Error 
bars are smaller than the symbol size. The inset shows the 
dependence on the spatial dimension of the exponent cp. At 
the upper-critical dimension of percolation, d = 6, the set of 
bridge bonds becomes dense having the spatial dimension. 



line. Above Pc^ as in the classical case, the percolation 
cluster is compact and there are no cutting bonds. For 
p < Pc, the set of cutting bonds is fractal with the same 
fractal dimension as the bridge-bond set, dcB = dsB, 
whereas at it is 1/i^. For the crossover scaling a similar 
ansatz to the one given by equation ([T]) is verified, where 
the argument of the scaling function is then {pc — p)L^ . 
The same hyperscaling of equation ([2| is obtained with 
if = d — dcB' At(i = 2, our numerical results corrobo- 
rate the hypothesis of the same C and (/p, for cutting and 
bridge bonds (see Appendix). For d > 2^ the set of cut- 
ting bonds is a line with dcB ^ and the one of bridge 
bonds has a dimension above d — 1, so that the fractal di- 
mensions differ. Above the critical dimension, the set of 
cutting bonds has dimension two, like the shortest path 
at Pc [32], for any p < Pc. 

Fisher [33| proposed a bond-site transformation to map 
bond percolation on a lattice onto site percolation on a 
covering graph. For example, the covering graph of the 
square lattice is obtained by adding all diagonal edges 
to every other face [34]. Considering such a mapping 
for ranked percolation and, following the analogy with a 
random landscape, where sites are sequentially removed 
from the lowest number (height) and suppressing global 
disconnection, the obtained cutting line is identical to the 
bridge one with the constraint applied in the perpendic- 
ular direction. The same is also observed for site-ranked 
percolation on a square lattice. Given the relation be- 
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FIG. 6. Following the rank, bonds are sequentially removed 
except the cutting bonds. Size dependence of the number of 
cutting bonds, Ncb^ for different fractions of occupied bonds 
p, namely p — pc — (circles), p — p — 0.49 (stars), and 
p = 0.2 (triangles). The solid lines are guides to the eye. At 
p — Pc^ NcB ^ L^/^, where v — 4/3 is the critical exponent 
related with the correlation length. For p < Pc, the number 
of cutting bonds scales with L^^^ . A crossover between the 
two regimes with the system size is observed (stars) for p 
in the neighborhood of pc. Systems of size have been 
considered, with L ranging from 8 to 512. Results have been 
averaged over 10^ samples for most system sizes and over 50 
samples for the largest one. Error bars are smaller than the 
symbol size. The inset shows the number of cutting bonds, 
NcB: as a function of occupied bonds, p, for two-dimensional 
lattices with L ranging from 64 to 256 (averaged over 10^ 
samples). The scaling is given by Eq. (1) in the article, where 
the argument is then {pc — p)L^ , with — 0.93. We obtain 
C = 0.56 =b 0.08, which is compatible with the one for bridge 
bonds. 



tween connectivity and disconnection, in both cases, the 
cutting version needs to be defined on different topolo- 
gies, namely, the star lattice for the square (on sites) and 
the same lattice for its covering graph, but with swapped 
cells of diagonal bonds. For the triangular lattice, the 
relation between cutting and bridges sites is straightfor- 
ward without requiring additional connections. 



II. DISCUSSION 

For several different models in 2D the same fractal di- 
mension as the one of bridge bonds has been reported. 
Here we discuss how some of them can be related with the 
process of fracturing ranked surfaces. Let us construct 
ranked percolation configurations on a random landscape 
starting with the bond with the largest number (height) 
and then occupying bonds sequentially in order of de- 
creasing number. Each time a chosen bond closes a loop, 
it is not occupied (loopless). One stops the procedure 



FIG. 7. Size dependence of the number of bridges, Nb^ on 
several different lattices. The solid lines are guides to the eye, 
all with slope 1.215. The best fit to each data set gives the 
same fractal dimension of 1.215 =b 0.005. Systems of linear 
size L have been considered, with L ranging from 32 to 4096. 
Results have been averaged over 10^ samples. 



when, for the first time, a path of connected sites spans 
from one side to the other (percolation threshold). Since 
all clusters are trees, the backbone and the shortest path 
are the same and their fractal dimension is identical to 
the bridge-bond line, as verified in Fig. [5] In fact, given 
the connectivity /disconnection relation discussed in the 
previous section, this line corresponds to the bridge line 
when bonds are occupied from the lowest to the largest 
number in the covering graph. This line also corresponds 
to the line of cutting bonds when bonds are removed 
from the smallest to the largest number. A similar pro- 
cedure was also proposed to obtain the minimum span- 
ning tree (MST), for which the same fractal dimension is 
found for the paths between all pairs of sites [19 . There 
bonds are sequentially removed (from the highest to the 
lowest value) under the constraint that all sites remain 
connected. If the constraint is relaxed such that only con- 
nectivity between two opposite borders is imposed, the 
cutting-bond line is obtained. Therefore, the cutting- 
bond line is the path between borders on the MST for 
which the largest value (height) is minimum, a valid re- 
lation in any spatial dimension d. 

The fractal dimension dsB was also observed for the 
backbone of the optimal path crack (OPC) [14^ 15 . For 
a random landscape, the sequence of optimal paths be- 
tween opposite borders is obtained and their highest site 
removed. Every such path crosses the bridge line. The 
highest site on successive paths is either on the bridge line 
itself or is higher than the lowest bridge-line site. There- 
fore, as in the loopless percolation described before, the 
removed sites percolate when all sites on the bridge line 
are removed, which is the backbone of the crack, giving 
the same fractal dimension dsB- This is in fact the case 



6 




FIG. 8. Number of bridge bonds, Nbb, as a function of the 
fraction of occupied bonds, p, for 3D (simple-cubic lattices) 
with different sizes L = {256, 512, 1024}. The scaling function 
proposed in this article is applied, with = 1.36, obtaining 
C = 1.0 =b 0.1. In the inset, Nbb has been rescaled by L^^^, 
where dsB = 2.50. All results have been averaged over 10^ 
samples. 



FIG. 9. Size of the largest cluster of sites on the edge of 
bridge bonds, Mbb, rescaled by L^^^, as a function of the 
fraction of occupied bonds p. The inset shows x defined by 
Eq. (Dl ). Results have been obtained on a square lattice with 
L ranging from 256 to 4096. All results have been averaged 
over 10^ samples. 



for the crack of any sequence of self-avoiding paths be- 
tween opposite borders. Since in this case the duality 
between cutting and bridges is not used (only valid in 
2D), the relation between OPC and bridge bonds is still 
true in higher dimensions. 

Let us now define the optimal minimax path (OMP) 
in the following way. One starts selecting the set of paths 
on the landscape for which the highest-ranked site is min- 
imum (minimax paths). Such set is then reduced to in- 
clude only the paths for which the second highest site is 
also minimum and one proceeds iteratively to the follow- 
ing sites, until a unique path is obtained. This path is 
the optimal path in strong disorder [16[|T7], since under 
such disorder strength each site is higher than the sum 
over the height of all sites with lower rank. This path 
is also identical to the backbone of the loopless-ranked- 
percolation cluster when occupied from the lowest to the 
largest number and, therefore, the cutting line in any 
dimension. 

In summary, suppressing connectivity (disconnection) 
between opposite borders on ranked surfaces leads to a 
fractal set of bridge (cutting) bonds, with a universal 
fractal dimension, even far away from the critical point 
of classical percolation. In 2D, there is an equivalence 
between bridges and cutting bonds and dsB = dcB^ 
whereas for d > 2 cutting bonds are still a fractal line 
but bridge bonds form a surface. The discussed mod- 
els are then split into two groups. The ones suppressing 
connectivity (e.g., watershed [lOl [12] and optimal path 
crack [HI [15]) are in the bridges universality class, while 
the ones keeping connectivity (e.g., optimal path on the 
MST or in strong disorder media fT6|) are in the univer- 



sality class of cutting bonds. For d > 6, dsB = d and 
dcB = 2, so we conjecture that the upper-critical dimen- 
sion of bridges and cutting bonds is also dc = Finally, 
we show that, at the percolation threshold of classical 
percolation, ranked percolation displays a theta-point- 
like crossover. 

This work opens up several challenges. Besides the 
need for a more precise numerical estimation of the 
bridge-growth exponent, it would be interesting to for- 
mulate a renormalization group scheme and to obtain the 
new set of exponents from analytic treatments such as, 
e.g., exact results in the mean-field limit and a Schramm- 
Loewner evolution in two dimensions. Another interest- 
ing possibility is to find the corresponding exponents in 
other related problems with different universality classes 
like, e.g., the Kasteleyn-Fortuin clusters in the g'-state 
Potts model [26l [35H38] with or without magnetic field. 
Regarding the fractal dimension of the surface of discon- 
tinuous percolation clusters [39, 40 , it would also be in- 
teresting to understand how it relates with the herein in- 
troduced bridge-bonds universality class. For the cutting 
bonds, the study of the crossover for higher dimensions 
represents another computational challenge. Finally, it 
would also be interesting to try to identify the third scal- 
ing field of our theta-like point. 



III. METHODS 

All numerical results have been obtained with Monte 
Carlo simulations. Results have been averaged over 10^ 
samples for 2D and 3D^ and 10^ samples for higher di- 
mensions. 
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Appendix C: Crossover scaling in 3D 

Figure [S] shows the crossover scaling for the number of 
bridge bonds Nb b on a simple-cubic lattice with different 
sizes. Bonds are occupied except for the bridge bonds. 
We applied the scaling function proposed in Eq. (1) of 
the article with = 1.36, giving ( = 1.0 ± 0.1. 



Appendix A: Cutting bonds 

In Fig. |6] we see the size dependence of the number 
of cutting bonds in 2D when bonds are removed with 
the constraint that cutting bonds are never removed. As 
discussed in the article for bridge bonds, a crossover is 
observed at p = between two fractal dimensions. In 
2D, for p < Pc the fractal dimension of the cutting bonds 
dcB is the same as observed for the bridge bonds for 
p > Pc^ whereas at Pc it is 



Appendix B: Universality of the fractal dimension of 
bridge bonds 

In Fig. [7| one sees results for the bridge-bond line in 
several different lattices, namely, square, star, triangu- 



Appendix D: Percolation of bridges 

As discussed in the article, the set of bridge bonds 
merges towards a single connected line. In Fig. [9] we 
see the size of the largest cluster of sites on the edge of 
bridge bonds Mbb^ rescaled by L^^^, as a function of p, 
for different system sizes. The inset shows defined as. 



1 



mIb 



(Dl) 



where the sum runs over all clusters of sites on the edge 
of bridge bonds and si is their size. Alike percolation on 
a line, in the thermodynamic limit, the bridge-bond line 
only percolates at p ^ 1. 
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